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ABSTRACT 


Two  series  of  experiments  were  performed  to  study  the 
interactions  between  both  the  wind-induced  drift  layer  and  the 
surface  wave  spectrum  and  a spatially  varying  subsurface  current 
field.  In  the  firsts  measurements  were  made  of  the  response  of 
the  drift  layer  and  surface  wave  spectrum  to  an  appreciable,  but 
unknown,  current  gradient.  This  was  accomplished  by  forcing  the 
current  from  a 90  cm-deep  diffuser  section  onto  a 68.7  cm-deep 
flat  beach.  The  flow,  thereby,  experienced  a current  gradient 
which  diminished  to,  and  remained  essentially  at,  zero  a short 
distance  downstream  of  the  beach  leading  edge.  r-JEije  results  in 
terms  of  the  velocity  difference  across  the  drift  layer  differ 
from  those  for  the  wind-only  case  and  will  not  acjmit  to  a simple 
Galilean  transformation  at  the  speed  of  th£-<f!Luid  below.  Pro- 
ceeding downwind,  however,  this  diffe^nce  decreases,  giving 
evidence  of  an  initial  phase  lag  iff  the  layer's  response  fol- 
lowed by  a relaxation  phenomenon.  ^Relaxation  times  inferred 
from  the  data  are  on  the  order  of  1 to  3 minutes.  For  .the  sec- 

ond  test  series,  the  beach  was  set  at  an  angle  of  2.65^  to  pro- 

“3  “*2  “1 

duce  strain  rates. iTT^he  range  10  < 3Uc/5x  < 10  (sec  ). 

S>Here  again,  the  results  indicate  that  the  response  of  the  drift 
layer  lags  behind  that  of  the  current,  as  strained  by  the  beach. 
The  measurements  of  surface  wave  spectra  taken  for  each  of  the 
conditions  tested  are  presently  being  analyzed  in  detail;  how- 
ever, some  preliminary  results  are  presented. 
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1.  INTRODUCTION 

Over  the  past  ten  years  or  so,  the  various  aspects  of  the 
problem  of  the  interaction  between  surface  waves  generated  under 
the  action  of  wind  and  subsurface  current  fields,  such  as  pro- 
duced by  a train  of  progressive  internal  gravity  waves,  have 
been  investigated  in  great  detail.  Areas  of  concentrated  effort 
have  included,  for  example,  the  generation  of  surface  waves  by 
wind,  the  dynamics  of  the  surface  waves  themselves,  the  mainte- 
nance of  surface  waves  by  wind,  the  dynamics  of  internal  waves, 
and  the  interactions  between  surface  and  internal  waves.  No 
doubt  major  steps  have  been  taken  both  experimentally  in  iso- 
lating the  phenomenon  and  studying  the  parts  played  by  the  vari- 
ous relevant  parameters,  and  theoretically,  in  modeling  the 
physics,  and,  thereby,  offering  predictions. 

Still,  the  broad  range  and  extremely  variable  nature  of  the 
flow  condition  encountered  in  the  oceans  make  the  problem  all 
the  more ' difficult  and  challenging.  Theoreticrl  models,  which 
describe  the  observed  interaction  between  wind-generated  surface 
waves  and  internal  waves  propagating  at  comparable  speeds,  fail 
when  the  speeds  are  strongly  disparate.  Several  mechanisms  have 
been  identified  and  each  has  its  own  merits;  however,  much  yet 
remains  to  be  done.  Some  highlights  of  the  currently  evolving 
models  of  and  pressing  questions  relating  to  wind-wave-current 
interactions  will  be  given  here  to  rationalize  the  motivation 
for  the  present  study. 
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Phillips^  first  considered  the  mechanism  of  blockage  to 
account  for  the  modifications  that  can  be  expected  when  a wind- 
generated wave  field  encounters  a surface  current  disturbance 
associated  with  a train  of  internal  waves.  He  showed  that  the 
surface  waves  will  suffer  significant  modifications  (or  block- 
age) when  the  propagation  speed  of  the  internal  wave  pattern  is 
approximately  equal  to  half  the  phase  velocity  of  the  incident 
surface  wave  train,  even  if  the  disturbance  velocity  in  the 
pattern  is  small.  This  effect  has  been  observed  experimentally 
in  the  laboratory  demonstrating  its  reality.  However,  the 
actual  flow  field  and  the  physics  of  wind-wave-current  inter- 
actions are  considerably  more  complex  than  admitted  by  Phillips ’ 
first  model.  In  particular,  no  account  was  taken  of  the  wind 
drift  layer  produced  by  the  shear  stress  that  exists  at  the 
air-water  interface. 

• 

When  the  wind  blows  over  the  ocean,  a shallow  boundary 
layer  flow  is  generated  at  the  water  surface.  The  effects  in- 
volving this  wind-induced  drift  layer  are  considered,  generally, 
to  be  of  two  types:  (1)  Those  in  which  the  drift  layer  changes 

the  surface  wave  dispersion,  and  (2)  Those  in  which  the  drift 

layer  itself  is  modified  by  the  current  field.  Modeling  the 

(2) 

former  case,  Phillips'  ' predicted  that  the  presence  of  the 
drift  layer  will  change  the  speed  of  propagation  of  even  those 
surface  waves  whose  wave  lengths  are  large  compared  to  the  layer 
thickness.  Changes  in  propagation  speed  in  turn  modify  the  in- 
teractions between  the  surface  waves  and  a perturbing  current 
field . 
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For  the  latter  case,  i.e.,  distortions  of  the  drift  layer 
itself,  modulations  of  the  surface  wave  field  may  be  produced 
indirectly  via  modification  of  the  drift  layer.  In  this  regard, 
Phillips^  has  derived  a quantitative  expression  for  the  drift 
layer  distortion  that  can  occur  when  the  length  scale  of  the 
drift  layer  is  small  compared  to  that  of  the  spartial  variation 
of  the  perturbing  current.  His  results,  in  terms  of  the  velocity 
jump  across  the  drift  layer  (i.e.,  the  surface  drift  relative  to 
the  fluid  below  the  layer),  indicate  that  severe  distortions, 
such  as  stagnation  points  on  the  surface  and  recirculating  re- 
gions below,  can  be  produced.  The  effects  of  viscosity,  both 
molecular  and  turbulent,  have  been  neglected  in  this  analysis. 
Consequently,  the  drift  layer’s  response  to  current  perturba- 
tions (on  a co-ordinate  system  moving  with  the  internal  wave 
pattern)  must  be  steady  in  time  and  always  in  phase. 

Vaglio-Laurin^  has  also  considered  the  distortions  pro- 
duced in  the  wind  drift  layer  due  to  a perturbing  current.  In- 
cluding the  effects  of  viscosity,  he  proposes  that  the  response 
of  the  drift  layer  to  fluctuations  of  the  current  will  exhibit 
phase  lags  and  relaxation  which  are  characteristic  of  diffusive 
processes.  Furthermore,  significant  distortion  to  the  drift 
layer  will  be  confined  within  a fraction  of  the  layer's  thick- 
ness below  the  surface.  Vaglio-Laurin' s model  draws  inspiration 

(5) 

from  an  earlier  work  by  Ligh thill  w on  a laminar,  viscous,  unsteady 

(6) 

boundary  layer  and  from  a study  by  Reynolds  and  Hussain  on  the 
mechanics  of  an  organized  wave  in  turbulent  shear  flow. 
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(3>) 

Recently,  Phillips'  ' has  developed  a theoretical 

formulation  for  the  modulations  produced  on  short  gravity  waves 

by  long  waves,  which  are  themselves  being  modulated  by  an 

internal -wave-induced  current  field.  The  problem  of  short  wave 

modulation  and  breaking  in  the  presence  of  long  waves  had  already 

(7) 

been  taker  up  by  Phillips  and  Banner'  ' . They  showed  that  non- 
linear modification  of  surface  drift  occurs  near  long  wave 
crests  and  troughs  to  modulate  the  drift  experienced  by  super- 
imposed short  waves.  One  very  important  consequence  of  this 
effect  would  be  to  reduce  the  maximum  amplitude  that  short  waves 
can  attain  when  they  are  at  the  point  of  incipient  breaking  near 
long  wave  crests.  Phillips  results  suggest  that  the  surface 
drift  over  a long  wave  can  vary  with  respect  to  the  phase  of  an 
attending  internal  wave  (because  the  long  wave  amplitude  varies) 
by  as  much  as  ten  percent  of  the  long  wave  speed.  While  only  a 
small  percentage  of  the  long  wave  speed,  this  modulation  repre- 
sents a substantial  fraction  of  the  phase  speed  of  the  short 
waves  riding  on  the  long  waves.  Thus,  significant  changes  of 
the  energy  density  and  propagation  speed  of  short  waves  may  be 
produced.  The  effects  of  viscosity  have  again  been  neglected 
so  that  perturbations  in  the  current  field  will  presumably  be 
felt  coherently  across  the  local  drift  layer. 

On  the  experimental  side,  the  study  by  Lewis,  Lake,  and 
Ko^  on  the  interaction  of  monochromatic  surface  and  internal 
waves  in  a two  fluid  system  is  noteworthy.  It  provides  valuable 
evidence  in  support  of  the  resonance  or  blockage  phenomenon  as 
previously  formulated  by  Phillips^^.  However,  no  wind  (and. 
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hence,  no  wind  generated  surface  drift  layer)  was  included.  The 
need  for  direct  experimental  measurements  of  wind-wave-current 
interactions  is  clear,  for  verifying  as  well  as  refining  the 
results  of  the  various  theoretical  models.  In  such  interactions, 
what  distortions,  if  any,  are  produced  in  the  wind  drift  layer? 

Is  the  process  diffusive  or  inviscid  in  nature?  How  is  the  sur- 
face wave  spectrum  modified  and/or  modulated  by  the  presence  of 
a subsurface  current  field?  These  questions,  among  those  raised 
in  a series  of  meetings  conducted  by  ARPA  last  spring  to  co- 
ordinate theoretical  and  experimental  efforts,  are  those  to  which 
the  experimental  test  program  reported  herein  has  been  addressed. 

The  experimental  test  program  was  conducted  in  two  over- 
lapping phases.  In  the  first  phase,  modulations  of  the  thin 
wind-drift  layer  by  a perturbing  current  were  measured  by  timing 
neutrally  buoyant  floats  of  several  sizes  between  two  fixed 
stations  in  the  wind  direction  under  the  various  test  conditions. 
In  the  second  phase,  the  effects  of  these  interactions  on  the 
surface  wave  slope-spatial  frequency  spectra  were  measured  util- 

(Q) 

izing  the  RRI  TV  camera  system^  ' . Simultaneously,  capacitance 
wave-height  probes  were  used  to  measure  surface  wave  height- 
temporal  frequency  spectra.  The  results  for  the  surface  wave 
slope-spatial  frequency  spectra  measurements  will  be  reported 
on  separately  by  Riverside  Research  Institute. 

Section  2 contains  a description  of  the  experimental 
facility,  the  measurements  taken  and  the  operating  characteris- 
tics. This  same  facility  has  been  utilized  in  related  studies 
over  the  past  several  years  with  only  minor  modifications. 
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There  exists,  therefore,  a large  amount  of  data  and  calibrations 
describing  the  operating  conditions  in  the  facility,  given  cer- 
tain set-up  parameters  such  as  fan  RPM,  fetch,  etc.  We  have 
drawn  upon  this  existing  information  whenever  possible,  encour- 
aged by  spot  checks  and  agreement  with  the  present  data. 

In  Section  3,  the  experiments  are  described  and  the  results 
are  presented.  The  important  question  of  experimental  accuracy 
for  drift  current  measurements  has  been  taken  up  as  a necessary 
part  of  the  present  study.  The  results  in  the  form  of  upper 
bounds  on  the  expected  standard  deviations  of  the  measurements 
are  also  presented  here. 


The  results  of  measurements  are  discussed  and  summarized  in 
Section  4.  Concluding  remarks  are  contained  in  Section  5,  along 
with  suggested  guidelines  for  further  investigation. 


2.  EXPERIMENTAL  FACILITIES  AND  MEASUREMENTS 


2 . 1 Wind-Wave -Current  Facility 


| 

' 

J 

[ 


The  HYDRONAUTICS  wind -wave-current  facility,  shown  schemat- 
ically in  Figure  1,  is  basically  a 1.5  m-wide,  1.55  m-deep,  and 
22  m-long  water  tank.  An  axial  flow  fan,  which  is  driven  at  an 
adjustable  speed  by  an  electric  motor  through  a magnetic  clutch, 
is  located  at  the 'upstream  end.  To  avoid  the  problem  of  wave 
reflections  that  would  occur  in  any  tank  of  finite  length,  a 
permeable  wave  absorber  has  been  installed  at  the  downstream 
end.  A removable,  sectional  cover  has  been  placed  on  the  tank 
to  create  a wind  tunnel  31  cms-high  over  a water  depth  of  about 
124  cms.  Wind  speeds  up  to  15  m/sec  can  be  generated  over  the 
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water  surface.  The  facility  is  also  equipped  with  a variable 
current  generation  system  also  shown  in  Figure  1.  The  system 
consists  of  a false  bottom  or  a submerged  beach  and  a recircu- 
lating pump  system,  including  a reversible  impeller  and  appro- 
priate ducting. 

The  primary  flow  variables  in  the  facility  are  the  wind 
speed  and  the  artificial  current  direction,  speed  and  streamwise 
gradient.  The  wind  speed  is  controlled  by  the  RPM  of  the  fan. 
The  current  direction  is,  of  course,  determined  by  the  sense  of 
rotation  of  the  impeller.  The  current  gradient  is  controlled 
by  the  inclination  of  the  beach  as  well  as  the  magnitude  of  the 
current,  while  the  current  is  directly  related  to  the  impeller 
RPM  and  the  inclination  of  the  beach.  It  is  the  current  grad- 
ient that  enables  the  simulation  of  certain  key  aspects  of  a 
moving  current  using  a stationary  current  system,  in  that  an 
equivalent  internal  wavelength  or  time-scale  for  a propagating 
current  can  be  defined.  It  should  be  noted,  however,  that  all 
measurements  are  hereby  made  relative  to  a coordinate  system 
that  effectively  "rides"  on  a stationary  internal  wave.  This 
point  must  be  considered  carefully  before  comparing  the  experi- 
mental laboratory  data  to  other  data  or  to  theoretical  predic- 
tions . 

2 . 2 Measurements 
2.2.1  Wind 

The  wind  velocity  profiles  in  the  tunnel  have  been 
previously  determined  by  vertical  traverses  of  a pitot-static 
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tube  supported  on  a precise  motor-driven  mechanism^10).  The  tube 
was  driven  through  the  air  stream  at  a constant  speed  of  2.5 
cm/sec  to  measure  the  wind  boundary  layer  near  the  water  sur- 
face. A differential  pressure  transducer  was  used  to  sense  the 
velocity  head.  Near,  but  not  overly  close  to,  the  water  surface 
the  wind  velocity  profile  was  found  to  follow  the  logarithmic 
law.  Friction  velocities ^and  roughness  lengths  were,  thereby, 
inferred  for  many  wind  speeds  and  are  presented  in  Figure  2. 

No  such  profiles  were  taken  in  the  present  experiments.  The 
presence  of  a small  (relative  to  the  wind  speed)  favorable  or 
adverse  current  no  doubt  effects  the  wind  profile,  and  the  value 
of  U*,  as  a modification  of  the  boundary  condition  at' the  air- 
water  interface.  For  small  currents,  this  effect  is  felt  to  be 
small  in  regard  to  the  establishment  of  U*  by  the  wind,  even  as 
the  wave -dependent  Stokes  drift  current  represents  only  a small 
percentage  of  the  wind-induced  drift  current  in  the  fetch- 
limited  labora^ory^"1) . In  other  words,  the  principal  effect  of 
a small  subsurface  current  is  felt  to  manifest  itself  as  a mod- 
ification of  the  surface  wave  spectrum  with  negligible  change  in 
the  shear  velocity.  As  the  magnitude  of  the  current  is  increased, 
however,  the  value  of  U*  increases  or  decreases  depending  on 
whether  the  current  is  adverse  or  favorable.  This  point  of  view 
is  supported  by  the  results  of  the  present  measurements  within 
the  drift  layer  and  will  be  taken  up  again  in  Section  4.  The 
results  of  Figure  2 were  accepted  as  being  approximately  appro- 
priate to  the  flow  conditions  of  the  present  experiments. 
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2.2.2  Surface  Wave  Height  and  Frequency 

A capacitance  probe  was  used  for  wave  measurements. 

The  probe  is  made  of  a 0.13  mm-diameter,  partially-submerged, 
Formvar- coated  (a  type  of  shellac  that  has  good  strength  and 
electrical  properties)  copper  wire  acting  as  one  electrode  and 
a fully- submerged  aluminum  plate  as  the  other.  The  change  in 
capacitance  of  the  wire  with  its  depth  of  submergence  is  con- 
verted through  an  AC  bridge  circuit  into  a change  in  voltage, 
which  can  be  maintained  in  a linear  proportion  over  a useful 
depth  range  with  suitable  electronics.  Within  the  frequency  and 
spatial  operating  range  of  the  probe,  the  output  represents  a 
trace  of  surface  height  as  a function  of  time  at  a fixed  point 
in  space.  It  is  possible,  then,  with  the  time  history  of  sur- 
face height  from  a single  probe  (which  in  reality  has  been  av- 
eraged over  two-dimensional  wave  number  space)  to  generate  a 
wave  height-temporal  frequency  spectrum.  The  probe  output 
voltage  was  fed  into  an  FM-tape  recorder  for  storage  to  be  pro- 
cessed at  a later  time  or  simultaneously  into  a visicorder  for 
an  immediate  graphical  display. 

Previous  studies  by  Wu^10^  of  the  surface  wave-height  and 
frequency  of  the  average  or  dominant  waves  under  various  wind 
conditions  did  not  include  the  effects  of  subsurface  current 
fields.  His  results  are  presented  here,  however,  in  Figure  3 
to  be  used  for  comparison  with  the  present  data  later  on  in  this 
report . 


mm 
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2.2.3  Surface  Drift  Current 

The  drift  current  immediate  to  the  water  surface  was 
measured  by  repeatedly  timing  floats  of  various  sizes  between 
two  stations  in  the  wind  direction  and  then  averaging  over  the 
results.  Spherical  particles  with  diameters  of  1.90,  3.18, 

6.35,  12.70  mm,  and  a specific  gravity  of  0.95  were  used  as 
floats.  The  velocity  of  each  float  was  taken  as  the  drift  cur- 
rent  at  the  depth  of  the  centroid  of  the  longitudinally  projec- 
ted area  of  the  submerged  portion  of  the  float.  For  these 
floats  used  in  the  experiment,  the  corresponding  centroid  depths 
have  been  estimated  (geometrically)  to  be  0.71,  1.19*  2.37, 

4.72  mm,  respectively.  In  some  cases,  two  additional ■ submerged 
floats  (each  constructed  of  a triangular  wooden  disk  and  a lor- 
mal  metal  wire  stem  of  the  correct  length  to  make  the  floats 

neutrally  buoyant  at  the  desired  depth,  see  Reference  8)  with 

• 

centroid  depths  of  12.7  and  24.5  mm  were  also  used.  For  many 
of  the  wind-current  set-ups  tested,  use  of  these  two  deeper 
floats  was  impractical  if  not  impossible.  In  particular,  for 
cases  with  adverse  subsurface  currents,  the  drift  layer  profile 
connecting  the  downwind  surface  drift  current  to  the  upwind 
subsurface  current  has  a zero  crossing,  typically  in  the  depth 
range  12.7  to  25.4  mm.  Under  conditions  of  near-zero  mean 
velocity  and/or  strong  velocity  gradients,  the  submerged  floats 
did  not  perform  well,  and,  therefore,  were  not  used. 

The  efficiency  of  the  timed-float  technique  has  been  dem- 
onstrated in  many  studies  of  the  drift  layer,  including  the 
present.  Albeit  cumbersome,  time  consuming,  and  unesthetic,  no 
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better  technique  has  been  discovered  as  yet.  Results  of  previous 
drift  current  measurements  by  Wu^^  for  the  case  without  subsur- 
face currents  are  presented  in  Figure  4;  included  here,  again, 
for  the  sake  of  comparison  to  the  present  data  taken  both  with 
and  without  currents.  The  important  questions  of  accuracy  and 
reoroducibility  of  results  obtained  with  the  timed-float  tech- 
nique will  be  considered  in  detail  in  Section  3.1.1. 

% 

2.2.4  Subsurface  Currents 

Subsurface  currents  produced  by  the  current  generation 
system  were  measured  by  timing  the  transport  of  nearly  vertical 
dye  streaks  between  two  stations.  Powdered  dyes*  Rhodamine  (red), 
Alphazurine  (blue),  and  Uranine  (green),  were  mixed  in  cold  tap 
water  and  simply  poured  into  the  stream  at  least  one  meter  up- 
stream of  the  measuring  section.  Observations  indicated  that 
below  the  first  2 to  4 cms  of  depth,  wherein  the  dye  dispersed 
very  quickly  owing  to  high  turbulence  levels,  and  extending  down 
to  at  least  30  cms,  the  time-mean  current  was  nearly  uniform. 
Whenever  possible,  measurements  were  also  taken  of  the  transit 
time  of  small  neutrally  buoyant  dirt  particles  in  the  flow  be- 
tween the  same  stations.  In  general,  the  results  were  in  good 
agreement.  As  with  the  drift  layer,  measurements  were  repeated 
several  times  and  'averaged  for  both  dye  and  dirt  particle  tech- 
niques . 

The  current  flow  below  the  surface  drift  layer  was  also 
quite  unsteady,  though  not  as  turbulent  as  the  surface  layer, 
and  under  some  of  the  conditions  tested,  exhibited  a three- 
dimensional  structure  in  the  form  of  longitudinal  and  vertical 
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vortex  tubes „ These  distortions  were  found  to  increase  with 
increasing  magnitude  of  current,  to  be  worse  for  adverse  currents 
(down  the  beach)  than  for  favorable  currents  (up  the  beach),  and 
to  decrease  with  distance  from  the  leading  edge  of  the  beach. 

All  of  these  symptoms  point  to  a poorly  diffused  and  straightened 
recirculating  current  flow.  In  an  investigation  aimed  at  mea- 
suring a particular  possible  source  of  flow  distortion,  extran- 
eous sources  should,  clearly,  be  avoided,  or  at  least  minimised. 

Thus,  the  recirculating  current  flow  was  considered  to  be  inad- 
equately diffused  and  straightened  in  the  present  set-up  to 
provide  reasonably  smooth  current  fields  for  currents  in  excess 
of  (approximately)  ±15  cm/sec. 

3.  THE  EXPERIMENT 

I 

3 . 1 Preliminary  Considerations 

I 

3.1.1  Accuracy  of  Drift  Layer  Measurements 

The  efficiency  of  the  timed-float  technique  for  mea- 

1 

suring  the  wind  drift  profile  has  been  demonstrated  in  several 
previous  HYDRO NAUTICS 1 reports  and  has  been  accepted  and  suc- 
cessfully employed  in  several  other  outside  experimental  facil- 
ities. However,  for  the  present  experiment  it  was  necessary  to 
measure  the  change  in  the  drift  with  and  without  a perturbing 
current  field.  To  meet  the  obvious  questions  of  accuracy  and 
statistical  confidence  required  for  meaningful  results,  system- 
atic experiments  were  first  performed  to  optimize  sample  size 
and  travel  distance  for  minimum  standard  deviations.  We  expect 
for  the  case  of  zero  current  gradient  that  increasing  either 
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sample  size  travel  distance  (for  floats  that  travel  in  a rea- 
sonably straight  path  between  the  measuring  stations)  will  de- 
crease the  measurement  standard  deviation.  With  a non-zero 
current  gradient,  on  the  other  hand,  the  travel  distance  has  to 
be  below  some  value  for  meaningful  estimates  oJ'  the  local  drift 
velocity.  For  such  tedious  and  time  consuming  measurements,  some 
attempt  at  optimization  was,  clearly,  in  order. 

To  accomplish  this  end,  float  measurements  for  a series  of 
relatively  stringent  flow  conditions  (with  zero  current  gradient) 
were  repeated  many  times  (as  many  as  120)  over  different  dis- 
tances of  float  travel  for  each  size  float  tested.  The  measure- 
ments were  then  grouped  into  different  sized  samples,  i.e., 
samples  per  reading,  and  averaged  to  give  a number  of  readings 
or  "measurements."  The  standard  deviation  of  the  "measurements" 
was  then  computed  and  compared. to  that  of  the  other  groupings. 

Put  in  more  concrete  terms,  the  question  we  sought  to  answer 
was  this  : Given  N samples  of  a measurement,  to  be  separated 

into  m groups  of  n samples  each,  and  averaged  over  n,  for  what 
value  of  n will  the  standard  deviation  of  the  m groups  be  mini- 
mum? This  study  could  more  properly  be  termed  an  "educated 
optimization"  in  that  "reasonable"  values  of  float  travel  dis- 
tance and  of  samples  per  measurement  (n)  were  selected,  based 
on  experience,  and  statistics  were  computed  and  compared  in  the 
vicinity  of  these  preselected  values. 

Trials  were  run  with  the  beach  in  the  horizontal  position, 
the  wind  speed  set  at  6.7  m/sec  (next  to  the  highest  speed 
tested),  and  for  subsurface  currents  of  0,  -1-9*2  cm/sec  (favorable) 
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and  -9.3  cm/sec  (adverse).  Float  diameters  (in  mm)  of  1.91, 

3 .18,  6.35,  and  12.7,  as  well  as  a 25.4  mm-deep  submerged  float, 
were  tested  over  travel  distances  of  .5,  .7,  and  1.2  meters. 

The  test  conditions  are  summarized  below: 

Beach  angle  = 0°  (zero  current  gradient) 

Wind  speed  = 6.67  m/sec 

Subsurface  current  = -9.3,  0.0,  +9.2  cm/sec 

Float  travel  distance  = .5,  .7,  1.2  m 

Float  diameter  = 1.91,  3.18,  6.35,  12.7  mm  plus 

25.4  mm-deep  float 
Trials,  ranging  from  15  to  .120  each. 

The  number  of  samples  per  measurement  (n)  was  taken  as  5,  10, 
and  15  (and  for  a few  cases,  20) . 

The  results  indicate  that  a float  travel  distance  of  1.2  m 
and  ten  samples  per  measurement  (the  high  and  low  are  thrown 
out  before  averaging)  yield  the  lowest  standard  deviations. 

These  are  presented,  in  part,  in  Table  1,  where  the  standard 
deviations  for  the  measurements  are  seen  to  vary  inversely  with 
float  diameter,  being  generally  larger  for  the  smallest  float. 
The  largest  standard  deviations  obtained  are  on  the  order  of  .5 
to  .6  cm/sec,  which,  relative  to  the  mean  velocity  estimates, 
corresponds  to  only  a 3 to  4 percent  variation.  It  should  be 
noted  that  the  statistics  no  doubt  vary  as  the  test  conditions 
are  changed.  However,  these  measurements  which  were  taken  with 
moderate  wind  speed  and  relatively  strong  currents  are  believed 
to  provide  reasonable  estimates  for  the  expected  measurement 
uncertainties  over  all  the  flow  conditions  tested,  including 
the  case  of  small  non-zero  current  gradient. 
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An  additional  factor  yet  to  be  mentioned  is  the  presmooth- 
ing or  filtering  built  right  into  the  raw  data.  During  the 
course  of  the  experiments  (in  fact,  for  all  experiments  using 
the  timed-float  technique),  those  particles  that  pop  out  of  and 
skip  along  the  water  surface,  visibly  dive  down  to  lower,  slower 
layers,  or  veer  from  the  wind  direction  along  the  center  of  the 
tank  by  more  than  ±15  cms  over  the  float  travel  distance,  are 
ignored.  The  first  two  effects  are  more  common  with  a water 
surface  having  large  amplitude  waves,  not  necessarily  sharp 
crested,  and  are  apparently  related  to  the  flow  fields  induced 
by  the  waves  themselves.  The  latter  effect  is  felt  to  be  caused 
by  three-dimensional  structure  within  either  the  air  stream  or 
the  subsurface  current  field.  For  example,  large-scale,  air- 
current  fluctuations  in  the  laboratory  caused  by  the  ventilation 
system  or  a suddenly  opened  or  closed  door  would  manifest  them- 
selves as  vorticity  in  the  wind  stream.  The  three-dimensional 
structure  in  the  current  field  has  been  already  discussed  briefly 
in  an  earlier  section.  The  testing  procedure  during  these  trial 
runs  was  the  same  as  that  used  during  the  experiments;  the  re- 
sults are  believed,  therefore,  to  carry  over. 

The  results  of  this  preliminary  study  are  important  for  the 
following  two  reasons:  They  prescribe  an  "optimized"  test  pro- 

cedure by  specifying  float  travel  distance  and  number  of  float- 
trials  per  measurement.  They  provide  estimates  of  the  accuracy 
to  be  expected  for  each  diameter  float.  Criteria  of  success,  of 
course,  lie  with  the  quality  of  experimental  data.  As  shall  be 
discussed  later  on,  when  considering  the  measured  data  in  detail. 
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the  estimated  accuracies  appear  to  be  quite  realistic,  if  not  a 
little  conservative. 


3.1.2  Selection  of  Test  Conditions 


The  following  theoretical  guidelines  as  to  the  most 
useful  and  interesting  test-flov/  conditions  were  supplied  to  us 
after  the  aforementioned  ARPA  organization  meetings. 


a.  Wind  speedy  < 8 m/sec,  to  which  was  added  the 
lower  limit  of  2.85  m/sec  below  which  the  surface 
activity  is  small  and  difficult  to  measure 
quantitatively . 


b.  Subsurface  current,,  U , both  adverse  and  favorable 


U 


- 4 5 U 


D 


o 


c max 


* + 4 


where  is  the  wind- induced  surface  drift  veloc- 

o 

ity  and  U is  the  maximum  current  at  the 

J c max 

downwind  end  of  the  sloping  beach. 
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Furthermore,  the  larger  values  were  felt  to  be  of 
special  interest  to  enhance  the  anticipated  per- 
turbations and,  thus,  to  facilitate  the  measure- 
ments. As  previously  discussed,  the  magnitude  of 

U was  constrained  by  consideration  of  flow 

c max  J 

geometry  so  that  in  some  of  the  cases  tested  these 
limits  were  exceeded. 


Strain  rate. 


c max 

from  the  range  of  subsurface  currents 


It  is  clear 
selected  in 
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b.  that  the  resonance  phenomenon  will  probably  not 
be  an' active  mechanism  here. 


3.1.3  Matrix  of  Flow  Conditions  Tested 


To  date,  two  series  of  experiments  have  been  performed. 
In  the  first,  measurements  were  made,  at  two  downwind  stations, 
of  the  drift  layer's  response  to  an  apprecialbe  but  unknown 
current  gradient.  This  was  accomplished  by  fixing  the  beach  in 
the  horizontal  position  at  a depth  of  68.7  cm.  The  current  flow 
in  passing  from  the  90  cm-deep  diffuser  section  onto  the  beach, 
experiences  a current  gradient  which  diminishes  to,  and  remained 
essentially  at,  zero  a short  distance  downstream  of  the  beach 
leading  edge.  The  flat  beach  case  is  seen  to  simulate  an  in- 
finitely long  internal  wave  and  the  drift  current,  as  measured 
along  the  beach,  should  then  reflect  the  layer's  response  to  an 
impulsive  disturbance.  Initially,  data  were  obtained  at  sta- 
tions 2 and  6;  138  and  38l  cms  from  the  leading  edge,  respect- 
ively, by  measuring  between  stations  1-3  and  5-7  (see  Figure 
1) . Data  in  this  set  is  coded  and  referred  to  as  B01 . Another 
set  of  data,  as  part  of  the  series  for  zero  beach  angle,  was 
taken  at  stations  3 (202  cms  from  leading  edge)  and  6,  in  a man- 
ner similar  to  B01;  this  set  is  coded  as  B02.  Surtace  wave 
spectra  (both  temporal-wave  height  and  spatial-wave  slope)  were 
also  measured  in  the  latter  set,  B02. 


For  the  second  test  series,  the  beach  was  set  at  an  inclin- 


ation from  the  horizontal  of  2.65°  to  produce  strain  rates. 


( i/u  Wdu  /dx)  < 8 X 10“f(cm-1).  Both  drift  layer  profiles 

v c max 7 v c ' 
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and  surface  wave  spectra  measurements  were  taken  at  station  6 
only.  This  test  series  is  coded  as  B4 . 

The  matrix  of  the  flow  conditions  actually  tested  for  each 
series  is  presented  in  Table  2 (a,  b,  and  c). 

3 . 2 Experimental  Results 

3.2.1  Flat  Beach  Series  - Test  Set  B01 

The  results  of  measurements  within  the  drift  layer  at 
an  upstream  and  a downstream  station  along  the  flat  beach  for  a 
moderate  wind  speed  and  both  a favorable  and  an  adverse  current 
are  shown  in  Figure  5 along  with  standard  deviation  bars.  Pre- 
sentation is  in  terms,  of  q,  the  velocity  relative  to  the  fluid 
below,  i.e.,  q = U_  - U . (IL,  is  the  measured  local  drift  vel- 
ocity  and  Uc  is  the  measured  local  current.)  If  the  drift  layer 
were  everywhere  in  local  equilibrium  with  the  current  field, 
the  current  would  merely  translate  the  layer  in  the 
Galilean  sense.  The  data  would  thereby  be  collapsible  onto  one 
curve  in  terms  of  q.  The  results  of  Figure  5 indicate,  however, 
that  this  is  not  the  case. 

It  should  be  noted  that  the  error  bars  given  here  were 
found  to  be  quite  realistic  in  that  data  runs  repeated  several 
days  apart,  while  attempting  to  reconcile  obvious  distortions 
in  these  plots,  consistently  yielded  results  that  fell  within 
the  bars.  Reproducibility  of  data,  therefore,  while  supporting 
the  earlier  estimates  of  measurement  error  also  pointed  out 
certain  unexpected  distortions  in  the  meanflow.  These  distor- 
tions are  believed  to  be  the  result  of  three-dimensional  flow 
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structure,  as  previously  mentioned,  caused  both  by  an  inadequately 
diffused  and  straightened  recirculating  current  flow  and  by  flow 
separation  from  the  beach-support  struts. 

After  fairing  smooth  curves  through  the  data  at  each  sta- 
tion, as  shown  in  Figure  5,  differences  were  taken,  at  a fixed 
depth,  between  the  results  with  and  without  current  and  are 
plotted  versus  distance  along  the  beach  in  Figure  6.  Here  the 
differences  are  seen  to  diminish  in  the  downstream  direction. 
Interpreting  the  differences  as  phase  lags  in  the  drift  layer's 
response  to  the  current's  history  coming  onto  the  flat  beach, 
the  observations  strongly  suggest  a relaxation  phenomenon.  Es- 
timates, of  a relaxation  time  constant,  t , were  inferred  using 
the  results  of  Figure  6 and  the  relationship 


(q  - qD) 

downstream 


(q  - qo) 

upstream 


e 


-t/T 


where 

L u 5 

avg 

Ax  = distance  between  stations  2 and  6,  and 

U = average  of  the  mean  velocities  at  stations  2 and 
avg 

6 at  the  given  depth. 

These  are  presented  in  Table  3(a)  for  each  measurement  depth. 

The  estimates  of  relaxation  time  are  based  on  rather  small 
differences  between  measurements  taken  at  only  two  stations. 
Moreover,  the  results  are  also  quite  sensitive  to  the  curves 
faired  through  the  data  points.  These  statements  may  be  more 
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easily  understood  from  a careful  examination  of  Figures  5 and  6. 
Note  in  Figure  5 that  .estimates  of  t are  based  on  the  separation 
between  the  curves  at  station  (5-7),  Figure  5(b),  as  compared  to 
the  separation  at  station  (1-3),  Figure  5(a).  The  quantitative 
details  are  presented  in  Table  3(a)  and  shown  as  a function  of 
downwind  position  in  Figure  6.  The  points  in  Figure  6 are  con- 
nected by  straight  lines  merely  to  indicate  the  trend  of  the 
data.  Our  assumption  is,  of  course,  that  the  data  points  lie 
on  exponential  curves.  The  expected  standard  deviation  for  each 
drift  velocity  measurement  has  already  been  given  as  roughly 
±|  cm/sec.  If  we  apply  this  error  range  to  the  measurements 
listed  in  Table  3(a),  where  the  differences  on  which  t is  based 
are  seen  to  vary  from  about  1.4  to  2.9  cm/sec,  the  computed  mag- 
nitude of  t is  formed  to  be  uncertain  by  a factor  ranging  from 
\ to  10.  Smoothness  and  reproducibility  of  the  actual  data, 
however,  promise  results  that  are  much  more  reliable.  That  this 
is  so  is  demonstrated  below,  especially  when  comparing  results 
for  different  data  runs.  Still,  we  must  view  the  calculated 
magnitudes  of  a with  these  point  in  mind. 

The  results  in  Table  3(a)  indicate  that  at  a mean  depth  of 
0.71  mm,  at  a wind  speed  of  6.67  m/sec,  time  constants  of  48 
and  21  seconds  are  appropriate  for  an  adverse  current  of  9.3 
cm/sec  and  for  a favorable  current  of  9*2  cm/sec,  respectively. 
The  magnitude  of  t is  apparently  constant  with  depth,  except 
for  the  4.7  mm  position.  At  this  depth  (and  below),  the  accur- 
acy of  the  velocity  differences,  and  especially  of  the  mean 
drift  velocity,  used  in  the  calculations  falls  off.  Little 
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confidence  is  placed  on  the  estimates  of  t at  the  4.7  mm  depth 
and  below.  More  will  be  said  on  this  topic  in  the  next  section 
wherein  the  data  from  the  more  complete  flat  beach  test  set  B02 
is  presented. 

3.2.2  Flat  Beach  Series  - Test  Set  B02 

The  second  set  of  experiments  performed  with  the  beach 
in  its  horizontal  position  (B02)  is  distinguished  from  the  first 
set  (B01)  on  two  counts 

a.  The  upstream  measuring  station  was  moved  down- 
stream from  station  2 for  B01  to  station  3 for 
B02 . This  change  was  made  principally  to  accom- 
modate the  RRI  TV  system.  However,  the  move  was 
also  justified  to  provide  a longer  flow  transi- 
tion length  from  the  beach  leading  edge  for  the 
wider  range  of  conditions  tested  in  B02. 

b.  A period  of  3 months  elapsed  between  B01  and  B02 , 
during  which  time  the  beach  was  adjusted  to  a 
2.65°  inclination  for  Test  B4  and  then  returned 
to  the  horizontal. 

As  previously  stated,  the  results  for  each  experiment  have  shown 
themselves  to  be  reproducible  within  the  estimated  accuracy  of 
the  data.  Whenever  possible  (i.e.,  whenever  the  experimental 
set-up  conditions  were  the  same),  therefore,  the  results  of  set 
B01  and  B02  are  compared  directly  to  one  another  later  on. 

The  results  of  measurements  at  both  upstream  and  downstream 
stations  within  the  surface  wind  drift  layer,  in  terms  of  the 


HYDRONAUTICS,  Incorporated 


-22- 


relative  velocity  q,  are  shown  in  Figures  7,  8,  9>  and  10,  each 
for  one  of  the  wind  speeds  tested.  The  same  effects  are  observed 
here  as  with  the  B01  data:  The  drift  layer  will  not  admit  to  a 

simple  Galilean  transformation  at  the  speed  of  the  current  below. 
The  layer  is  apparently  lagging  the  current  field  and  relaxing 
to  steady  state.  The  flow  combinations  of  a 2.75  m/sec  wind 
speed  and  5.1  cm/sec  favorable  current,  shown  in  Figure  7,  and 
of  a 7.65  m/sec  wind  speed  and  a 22.5  cm/sec  favorable  current, 
shown  in  Figure  10,  yielded  highly  distorted  drift  layer  pro- 
files. All  of  the  other  conditions  tested,  however,  yielded 
profiles  of  the  now-f amiliar  drift  layer  form;  that  is,  mono- 
tonically  diminishing  to  zero  with  depth.  The  relative  velocity 
q is,  of  course,  dependent  upon  measurements  of  the  imposed  sub- 
surface current  field. 

Figure  11  shows  the  distribution  of  subsurface  current 
along  the  beach  as  determined  by  detailed  measurements  in  the 
centerline  plane  for  each  flow  condition  tested.  Here  U is 
the  local  value  of  the  subsurface  current,  which  was  nearly 
uniform  within  4 to  30  cms  of  depth  at  a fixed  streamwise  sta- 
tion. In  other  words,  Uc  is  the  hypothetical  1-D  current  along 
the  beach.  U is  the  arithematic  average  of  U measured  over 
the  length  of  the  beach.  It  should  be  noted  that  the  presence 
of  a wind  blowing  over  the  water  surface  had  a negligible  direct 
effect  on  the  subsurface  current  field.  The  small  variation 
with  streamwise  position  shown  in  Figure  11  is  believed  to  be 
caused  principally  by  3-D  structure  in  the  actual  recirculating 
flow. 
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Proceeding  along  the  lines  described  above  for  the  B01 
data,  smooth  curves  were  faired  through  the  data  points  (except 
for  the  two  peculiar  cases  mentioned  above),  differences  formed 
and  relaxation  times  estimated.  Table  3(b)  presents  the  results 
of  these  calculations.  Here  we  find  evidence,  as  in  Table  3(a) 
for  B01 , that  t is  constant  with  depth.  Here,  again,  the  re- 
sult s for  7.7  mm  depth  (and  below)  must  be  viewed  with  an  eye 
towards  the  previously  estimated  factor  of  uncertainty  for  r of 
from  § to  10.  We  expect,  on  physical  grounds,  that  the  time 
constant,  describing  the  relaxation  of  the  upper  layers  to  the 
perturbing  current  below,  should  decrease  with  depth.  The  re- 
sults in  Table  3(b)  further  suggest  that  the  magnitude  of  t 
increases  with  decreasing  wind  speed  and  decreases  with  favor- 
able current.  At  a mean  depth  of  0.7l  mm,  a wind  speed  of  7.75 

m/sec,  and  for  an  adverse  current  of  approximately  8.2  cm/sec, 

• 

the  magnitude  of  r is  on  the  order  of  200  seconds. 

Comparing  the  results  of  test  B01  to  B02  at  a wind  speed  of 
6.67  m/sec  and  an  adverse  current  of  (approximately)  9-0  cm/sec, 
the  latter  are  seen  to  be  larger  than  the  former  by  a factor  of 
2 to  3.  This  variation  is,  however,  well  within  the  range  of 
the  expected  uncertainty  of  a and,  in  fact,  provides  a valuable 
check  on  the  accuracy  of  the  results  between  separate  but, 
otherwise,  identical  experiments.  Alternately,  that  data  from 
the  two  experiments  (with  identical  flow  set-up  condition)  ct  Id 
be  taken  together  to  compute  a single  time  constant  based  on  a 
larger  sample  size.  An  examination  of  the  velocity  differences 
for  these  two  tests,  given  in  Table  3(a)  and  (b),  reveals  that 


a*. 


— ...  — 
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this  procedure  would  introduce  the  full  estimated  range  of  un- 
certainty to  the  results,  .in.  place  . .of  the  smaller- factor  when 
treating  the  tests  separately.  The  velocity  differences  are 
seen  in  Table  3 to  be  in  agreement  with  one  another  within  the 
expected  accuracy  of  ±-|  em/sec,  which  leads  to  the  stated  range 
of  uncertainty  for  t . This  result  is  not  so  surprising.  It 

points  out  that  data  taken  continuously  during  one  experiment 

% 

shows  slightly  less  variance  than  data  taken  over  separate  ex- 
periments (realizations).  The  stated  accuracies  of  velocity 
measurements  are  based  on  readings  taken  over  many  realizations 
and  are  accepted  as  the  expected  accuracy.  For  the  experiments 
in  question,  then,  there  is  nothing  to  be  gained  by  treating 
the  results  together  over  treating  each  data  set  separately  for 
the  time  constant  estimates. 

Wave-height  temporal  frequency  spectra  were  generated  from 
the  calibrated  output  signal  of  a capacitance  wave-height  probe 
utilizing  a SAICOR  (Model  SAI-52B)  Real  Time  Spectrum  Analyzer/ 
Digital  Integrator.  A five -minute  long  signal  record  could  be 
processed  by  utilizing  the  SAICOR' s integrator  section.  Thirty- 
two  separate  spectra  are  computed,  one  after  the  other,  while 
inputting  a five-minute  long  signal  record,  and  frequency  aver- 
aged by  integration  to  yield  a smoothed  output  spectrum  for  the 
full  record.  The  results  for  only  one  wind  speed  will  be  pre- 
sented at  this  time,  as  being  typical,  to  bring  to  light  some 
important  observations  in  the  data.  The  remainder  of  the  spec- 
tral data  and  a complete  analysis  will  appear  under  separate 
cover  in  the  near  future  to  report  the  final  results  ancl 
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interpretations  of  the  joint  HYDRONAUTICS,  Incorporated-Riverside 
Research  Institute  1974  tank  tests. 

Figures  12  and  13  show  spectra  measured  at  stations  3 and 
6,  respectively,  at  a constant  wind  speed  of  *1.75  m/sec.  The 
vertical  scale  here  is  decibels  for  convenience  of  interpreta- 
tion. Mote  that  10  decibels  is  equal  to  1 decade  for  power. 

For  clarity  of  presentation,  the  spectral  plots  tagged  as  2 and 

*'v 

3 are  drawn  with  their  frequency  axis  shifted  horizontally  rel- 
ative to  1 as  indicated  on  the  figures.  Comparing  the  results 
from  each  figure  for  the  wind-only  case,  the  high-frequency 
portion  of  the  spectrum  centered  on  the  first  harmonic  of  the 
dominant  wave  is  seen  to  attenuate  with  fetch.  This  is  an  ex- 
pected result.  At  both  stations,  the  effect  of  a favorable 
current  field,  compared  to  the  wind-only  case,  is  to  increase 
the  frequency  of  the  dominant  wave  and  to  significantly  decrease 
the  energy  (or  amplitude)  level  of  surface  wave  activity.  An 
adverse  current  field,  on  the  other  hand,  decreases  the  frequency 
of  the  dominant  wave  and  only  slightly  increases  the  energy 
level.  These  effects  will  be  further  discussed  in  Section  4. 

3.2.3  Inclined  Beach  Series  - Test  B4 

Because  of  considerations  of  time  and  completeness, 
measurements  were  taken  only  at  station  6 for  this  beach  geometry. 
We  felt  that  it  would  be  most  useful  and  expedient  to  put  our 
time  and  energy  for  these  first  tests  Into  a broad  range  of 
flow  conditions  at  one  station  to  isolate  the  discernable  ef- 
fects of  wind-wave -current  interactions  and  to  consider  carefully 
our  ability  to  quantitatively  measure  them.  Furthermore,  the 
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RRI-TV  system,  which  was  plagued  with  time  consuming  electronic 
and  mlcrophonic  problems  that  weighed  heavily  on  the  program 
time  budget,  is  not  eadily  movable. 

Figures  Hi,  15,  16,  and  17,  each  one  representing  one  of 
the  wind  speeds  tested,  present  the  results  of  measurements 
within  the  drift  layer.  As  for  the  flat  beach  data,  the  mea- 
surements are  expressed  as  relative  velocity  1 = UD  ‘ uc  > 
where  U is  the  local  current  speed.  The  same  general  trends 
are  evident  here  as  compared  to  the  flat  beach  case.  That  is, 
the  (relative)  adverse  currents  are  larger  than,  and  the 
(relative)  favorable  currents  are  smaller  than  the  wind-only 
currents ; again,  suggesting  a phase  lag  in  the  response  of  the 
drift  layer  to  perturbations  in  the  subsurface  current  field. 

Smooth  curves  have  been  faired  through  the  data  points  for 
the  two  lower  wind  speeds,  shown  in  Figures  14  and  15.  The  re- 
sults for  the  higher  wind  speeds  (Figures  16  and  If)  and  favor 
able  currents  are  not  so  easily  described  in  that  the  data 
points  are  scattered  and  cannot  be  connected  by  a simple  curve, 
in  all  cases  the  data  were  reproducible.  Therefore,  the  quali- 
tative differences  between  the  results  at  different  wind  speeds 
are  viewed  as  real  distortions  in  the  mean  flow  field  that  are 
time-independent.  Measurements  of  the  local  subsurface  current 
speed  were  also  taken  for  each  test  set-up  in  order  to  construct 

the  relative  velocity  q. 

Figure  18  is  a plot  of  the  experimental  subsurface  current 
speeds  measured  along  the  length  of  the  beach  in  the  longitud 
inal  plane  of  the  centerline.  Comparing  the  measurements  to 
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the  theoretical  1-D  result  (representing  a hyperbolic  variation 
of  current  with  local  water  depth)  in  Figure  18,  it  .is  clear 
that  the  experimental  f.l.ow  was  rather  more  complex.  The  effect 
of  three-dimensionality  of  the  current  field  on  the  drift  mea- 
surements is  unknown  at  this  time.  As  previously  discussed, 
the  issue  was  avoided  as  far  as  possible  by  limiting  the  range 
of  current  speed  to  minimize  3-D  distortions. 

Wave  height -temporal  frequency  spectra,  for  a wind  speed 
of  Jl.75  m/sec  in  combination  with  favorable,  adverse,  and  zero 
current  fields  on  the  inclined  beach,  are  presented  in  Figure 
19.  Here,  again,  the  spectral  plots  have  been  horizontally 
shifted  relative  to  one  another  for  clarity.  Comparing  the 
results  for  the  wind-only  case  in  Figure  19  (for  an  inclined 
beach)  to  those  in  Figure  13  (for  a flat  beach)  at  the  same 
wind  speed  and  fetch,  we  find  agreement.  Increasing  favorable 
current  results  in  a higher  dominant  wave  frequency  and  in-  a 
strikingly  lower  spectral  energy  level.  Increasing  adverse 
current  yields  a lower  dominant  wave  frequency  and  only  a very 
slightly  higher  spectral  energy  level.  Another  effect,  that  is 
not  so  obvious  in  Figure  19  in  which  all  three  spectra  are  plot- 
ted on  the  same  scale  to  highlight  relative  differences  in  the 
high  energy  portions,  is  that  an  adverse  current  increases,  and 
a favorable  current  decreases,  the  relative  wave  energy  in  the 
high  frequency  portion  of  the  spectrum  in  the  vicinity  of  the 
first  harmonic  of  the  dominant  wave.  This  point  is  illustrated 
in  Figure  20.  The  spectrum  for  the  adverse  current  run  (curve 
1,  Figure  19)  and  for  the  favorable  current  run  (curve  3,  Figure 
19)  are  scaled  to  yield  the  best  comparison  of  relative  energy 
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levels  and  plotted  together  in  Figure  20.  The  adverse  current 
run  is  seen  to  have  a dominant  temporal  frequency  of  3.08  cps 
and  a clearly  discernible  first  harmonic.  The  favorable  current 
run,  on  the  other  hand,  has  a dominant  frequency  of  ^1 . 51  cps  and 
a barely  perceptible  first  harmonic.  The  fetch  is  effectively 
reduced  by  the  favorable  current  and  increased  by  the  adverse 
current,  compared  to  the  wind-only  case.  In  section  3.2.2,  the 
effect  of  fetch  was  shown  to  reduce  the  energy  of  the  harmonic. 
Therefore,  the  spectrual  modification  noted  above  is  felt  to  be 
due  to  current  rather  than  to  fetch. 

i|  . DISCUSSION 

There  are  three  important  aspects  of  the  experiment  to  be 
further  discussed  in  this  section;  namely,  wind-wave  generation, 
mean  recirculating  flow-establishment  and  characteristics,  mea- 
surements of  wind-drift  layer  distortions  and  surface  wave 
spectra  modulations.  The  latter  is,  of  course,  the  most  rele- 
vant here.  However,  to  establish  the  highest  level  of  confi- 
dence in  the  experimental  results,  our  aim  is  to  describe  the 
test  conditions  and  to  present  the  details  of  the  measurements 
as  completely  and  as  accurately  as  possible. 

*1 . 1 Wind -Wave  Generation 

The  wind-wave  tank  has  been  utilized  over  the  past  several 
years  in  studies  directly  related  to  the  present  one.  We  have 
attempted  to  draw  from  the  large  body  of  calibration  data  and 
operating  characteristics  in  existence  for  the  facility  whenever 
appropriate  and  possible.  Previous  studies  differ  from  the 
present  one  in  at ‘least  two  important  respects.  In  the  present 
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case,  measurements  were  taken  at  several  downwind  stations, 
while  in  the  past,  attention  was  normally  concentrated  at  one 
station  or  at  least  over  only  a small  downwind  interval.  The 
present  results  include  the  effects  of  wind  and  of  current, 
while  previous  studies  generally  include  that  of  the  wind  alone. 
Figure  3 presents  the  'wave-characteristics  of  the  water  surface 
as  a function  of  wind  speed-  for  the  wind-only  case,  for  both  the 

'v 

previous  and  the  present  studies.  Similarly,  measurements  of 
the  wind-induced  surface  drift  current  have  been  presented  in 
Figure  H . The  important  points  to  note  here  are  that  while  the 
agreement  between  the  previous  and  the  present  data  are  fairly 
good  there  is  a discernible,  but  small,  fetch  effect  due  to  the 
downwind  evolution  of  the  surface  wave  spectrum  and  the  wind 
profile.  We  treated  the  results  as  being  essentially  indepen- 
dent of  this  effect  by  comparing  (at  each  station)  the  data  for 

• 

the  cases  with  non-zero  current  to  that  for  the  wind-only.  In 
other  words,  we  treated  the  wind-only  data  at  each  station  as  a 
"tare"  measurement  and  were  able  thereby  to  proceed  further  with 
the  data  to  estimate  the  relaxation  time  constants  as  previously 
described. 

The  results  for  the  two  lower  wind  speeds  and  corresponding 
current  fields  hav.e  been  given  for  the  flow  over  a flat  beach  in 
Figures  7 and  8 and  for  the  flow  over  an  inclined  beach  in  Fig- 
ures It  and  15.  Examination  of  the  data  in  each  of  the  figures, 
comparing  the  wind-only  case  to  that  for  both  favorable  and  for 
adverse  current,  it  appears  that  the  slopes  dq/dz,  in  the  linear 
region  near  to  the  water  surface,  are  the  same.  In  other  words. 
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indications  are  that  for  the  lower  values  of  wind  speed  andx 
currents  tests,  the  friction  velocity,  w*  = vw  (dq^/dz)  , 

where  vw  is  the  kinematic  viscosity  of  the  water,  is  "a  direct' 
function  of  wind  speed  only  and,  thereby,  also  a weak  function 
of  fetch.  This  is  not  so  for  the  higher  wind  speeds  and  cur- 
rents tested,  as  indicated  in  Figures  3,  9,  and  10  for  the  flat 
beach,  and  Figures  1.6  and  17  for  the  inclined  beach.  For  these 
cases,  the  slopes  of  the  drift  profiles  at  a given  wind  speed 
apparently  vary  with  the  imposed  current  field;  We  expect,  on 
kinematic  grounds  at  least,  that  the  shear  velocity  will  be  in- 
creased by  the  presence  of  an  adverse  current  field  and  decreased 
by  a favorable  one.  These  results  sup^gest  that  the  modification 
of  friction  velocity  by  a subsurface  velocity  is  a second-order 
effect  of  the  ratio  of  wind-induced  surface  drift  to  current 
speed  and  remains  small  except  as  the  ratio  becomes  unity  or 
larger.  Data  are  not  yet  available  to  resolve  this  point;  di- 
rect measurement  of  the  shear  stresses  may  be  necessary. 

In  summary,  then,  the  wind  and  surface  wave  characteristics 
measured  in  this  study  agree  fairly  well  with  those  of  previous 
studies.  There  is  a discernible  effect  of  fetch  on  the  results 
at  different  measuring  stations.  The  shear  velocity  is  seen  to 
be  decreasing  downwind . Taking  the  surface  drift  currents 
measured  at  upstream  and  downstream  stations  to  Figures  2 and 
>\  , we  can  estimate,  for  the  wind-only  case,  the  local  wind 
speed  and,  thereby,  the  magnitude  of  the  shear  velocity  at  each 
station.  Between  stations  2 and  6,  the  changes  in  shear  veloc- 
ity Inferred  in  this  manner  are  on  the  order  of  10  percent/meter. 
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The  presence  of  a subsurface  current  field  apparently  has  an 
effect  on  the  magnitude  of  the  shear  velocity  which  becomes  more 
pronounced  as  the  ratio  of  wind-induced  surface  drift  to  current 
speed  increases.  This  effect  is,  however,  not  yet  understood. 

l\ . 2 Mean  Recirculating  Flow  - Establishment  and  Characteristics 

The  submerged,  stationary  beach-recirculat ing  current  sys- 
tem is  an  essential  part  of  these  experiments  in  that  it  provides 
subsurface  currents  and  current  gradients  to  simulate  important 
features  of  progressive  internal  waves  in  a laboratory  facility 
of  limited  size.  The  present  configuration  has  been  shown  to 
be  incapable  of  producing  uniform,  mean  current  fields  over  a 
large  speed  range.  This  problem  can  and  will  soon,  no  doubt,  be 
resolved  with  the  addition  to  the  system  of  properly  designed 
diffusing  and  smoothing  sections.  In  order  that  this  source  of 
unwanted  flow  distortions  may  not  become  a rug  under  which  all 
difficulties  (including,  perhaps,  some  subtle  aspects  of  the 
physics)  and  unexplainable  data  be  swept,  it  is  appropriate  that 
we  discuss  it  further  with  all  the  facts  in  hand. 

Meandering  of  the  floats  from  their  drop-points  along  the 
tank  centerline  to  cither  side  wall  gave  the  first  indication  of 
non-parallel  current  flow.  Float  size,  beach  orientation,  fetch, 
as  well  as  wind  speed  and  current  were  all  important  parameters. 
Sometimes  a float  would  travel  along  the  centerline  for  a meter 
or  more  and  then  suddenly  veer  to  the  side.  Once  a float  left 
the  centerline  and  moved  to  one  side  or  the  other,  it  remained 
there  as  it  was  carried  downwind.  We  investigated  the  structure 
of  the  current  field  visually  by  means  of  vegetable  dyes  and 
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were  thus  able  to  discern  the  presence  of  longitudinal  vorticity 
within  the  current  field.  The  strength  of  this  vorticity,  as 
indicated  by  the  rotational  velocity  of  the  dyed  fluid  and  the 
coherence  of  its  3-D  structure,  varied  with  the  speed  of  the 
current  ('i.e.,  impeller  RPM)  and  the  orientation  of  the  beach 
to  the  current.  Based  on  dye  observations  and  consideration  of 
the  vorticity  theorems  of  basic  fluid  mechanics,  we  feel  that 
the  vorticity  we  saw  is,  for  the  most  part,  the  natural  outcome 
of  the  geometry  of  the  recirculating  flow  over  the  beach  and  of 
the  impeller  which  drives  the  flow.  Both  of  these  effects  can, 
of  course,  be  greatly  attenuated  if  not  eliminated,  with  proper 
smoothing  and  straightening  of  the  flow. 

An  additional  source  of  vorticity  and  three-dimensionality, 
not  yet  discussed,  results  from  the  flow  geometry  for  the  case 

of  adverse  currents  in  the  wtnd-wave-current  facility.  Referring 

* 

to  Figure  21,  the  drift  profile  for  the  case  of  adverse  currents 
must  have  a zero,  crossing  to  match  the  wind-induced  drift  layer 
to  the  current  field  below . Stagnation  points  would,  then, 
exist  at  the  surface  somewhere  both  upstream  and  downstream. 

It  is  obviously  impossible  to  maintain  the  "bubble"  shown  and 
satisfy  a two-dimensional  mass  balance.  The  flow  field  must  be 
three-dimensional  someplace,  at  least  near  the  stagnation  points. 
This  three-dimensionality  is  a form  of  vorticity. 

In  general,  distortions  to  the  mean  flow  were  found  to  be 
small,  and  the  occurrence  of  veering  floats,  infrequent,  for 
both  adverse  and  favorable  currents  ^.15  cm/sec.  Put  in  more 
quantitative  terms;  for  currents  £ 15  cm/see,  it  was  difficult 
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to  perceive  any  3-D  vertical  structure  in  the  dyed  (turbulent) 
flow  fields  and  we  could  expect  on  the  average  of  one  float  in 
10  to  veer.  By  comparison;  for  currents  in  excess  of  20  cm/sec 
we  could  expect  as  many  as  5 out  of  10  floats  to  veer.  The  dis 
tortlonS;  shown  in  Figures  10;  l6 , and  17,  are  believed  to  be 
the  result  of  the  presence  of  strong  longitudinal  vorticity  in 
the  current  field. 

Flow  separation  from  the  beach  support  struts  at  the  side 
walls  of  the  tank  were  also  troublesome.  These  were;  however; 
more  of  a problem  at  the 'lower  current  speeds  (and  correspond- 
ing wind  speeds);  first  of  all;  because  the  surface  wave  activ- 
ity and  the  turbulence  level  in  the  water  were  less,  enabling  a 
longer  life;  and  second  of  all;  because  they  could  more  readily 
move  away  from  the  wall  toward  the  center  of  the  tank  where  the 
drift  measurements  were  taken.  The  distorted  drift  profile; 
shown  in  Figure  7 for  a favorable  current  of  5.1  cm/seC;  is 
believed  to  be  the  result  of  vertical  vorticity  shed  from  the 
beach  support  struts.  The  solution  to  this  problem  would  be 
simply  to  redesign  the  beach  and  support  it  from  below. 

. 3 Measurements 

b .3.1  Wind  Drift  Layer  Distortions 

The  results  of  our  measurements  indicate  that  changes 
in  the  structure  of  the  wind-induced  drift  layer  produced  by 
the  simulated  Internal  wave  current  field  in  these  experiments 
are  discernible  and  measurable;  though  quantitatively  quite 
small.  It  is  the  interpretation  of  such  measurements  that  make 
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thern difficult  and  not  the  actual  data  taking  which  is  an  ob- 
jective process.  As  previously  discussed,  the  preliminary  ex- 
periments to  optimize  the  data-taking  scheme  for  minimum  stand- 
ard deviation  have  proven  to  be  very  helpful  in  this  regard. 

The  data,  bars,  shown  in  Figure  5,  that  are  believed  to  be 
appropriate  for  all  the  measurements,  are  the  outcome  of  this 
preliminary  study.  Moreover,  the  smoothness  and  self-consistency 
of  the  data--w.ith  exception  of  those  distortions  particular  to 
the  background  mean  flow  as  discussed  above--tend  to  increase 
our  confidence  in  the  results  as  tabulated  in  Tables  3(a)  and 
3(b).  The  same  qualitative  results  were  obtained  for  all  flow 
conditions  tested  and  for  both  the  flat  and  the  inclined  beach 
configurations;  that  is,  the  relative  velocity  in  the  presence 
of  an  adverse  current  field  was  always  greater  than,  and  in 
the  presence  of  a favorable  current  field  was  always  smaller 
than  that  for  the  wind  alone.  These  differences  apparently 
decrease  in  the  downwind  direction.  In  other  words,  some  time 
is  required  for  the  drift  layer  to  adjust  itself  to  changes  in 
the  velocity  of  the  fluid  below. 

An  explanation  for  the  smallness  of  the  measured  changes 

in  the  drift  layer  can  be  seen  from  the  results  in  Tables  3(a) 

and  3(b).  The  actual  travel  time  between  measuring  stations 

at  a given  mean  depth  below  the  surface  is  estimated  by 

Ax:/U  . Comparing  these  values  to  those  for  the  estimated 
avg 

relaxation  time  we  see  that  travel  times  were  typically 
only  10  to  20  percent  of  the  relaxation  times.  This  is,  of 
course,  another  way  of  saying  that  the  distance  between  the 
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measuring  stations  was  small  compared  to  the  relaxation  length 
scale.  These  results  also  suggest  that  for  a constant  wind 
speed,  the  ratio  of  a for  an  adverse  current  of  (approximately) 

9 cm/sec  to  t for  a favorable  current  of  the  same  magnitude  is 
(approximately)  2.  However,  this  ratio  for  an  adverse  current 
of  9 cm/sec  and  a favorable  current  of  23  cm/sec  jumps  to 
( app  roximately ) 10  . 

There  is  insufficient  information  at  this  time  to  permit 
parameterization  of  the  flow  processes.  Still,  certain  trends 
pointed  out  in  Section  3.2  are  evident  and  worth  mentioning 
here  again:  Relaxation  times  are  seen  to  be  nearly  constant 

with  mean  depth,  to  decrease  with  wind  speed,  to  - \crcase  with 
adverse  current  speed,  and  to  decrease  with  favorable  current 
speed.  Also,  the  effects  of  current  gradient  cannot  yet  be 
assess*,  d and  await  further  tests. 

4.3-2  Surface  Wave  Spectra  Modulations 

Awaiting  the  RRI  report  on  their  measurements  of  the 
wave  slope  spectra  for  the  flow  conditions  reported  here,  only 
preliminary  remarks  can  be  made  at  this  time  based  on  our  wave 
height  data.  Figures  12  and  13  for  the  flat  beach  tests  (B02) 
and  Figures  19  and  20  for  the  inclined  beach  tests  (B4)  have 
already  been  presented  as  typical.  A line  drawn,  on  the  figures 
at  a slope  of  -5  provides  for  comparison  of  the  high  frequency 
portion  of  the  spectrum  in  the  wave-wind  facility  to  that 
predicted  by  Phillips (12 ) for  the  case  of  infinite  fetch.  A 
recent  calibration  of  wave-height  probe  frequency  response 
shows  a fall-off  of  -20Db  per  decade  from  a frequency  of  4 to  8 
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cps.  Correcting  the  spectra  already  presented  for  this  effect 
results  in  fairly  good  agreement  with  Phillips'  (-5)  prediction 
in  the  high  frequency  portion  of  the  spectrum.  Corrected  spectra 
are  not  shown  here  hut  will  be  presented  in  detail  in  the  afore- 
mentioned final  report  of  the  HYDRONAUTICS,  Incorporated- 
Riverside  Research  Institute  1974  tests. 

The  spectra  measured  in  the  limited  fetch,  wind-wave  tank 
relative  to  Phillips'  (-5)  spectrum  are  over-driven  at  the  fre- 
quency of  the  dominant  wave,  exhibit  a second  peak  in  the 
vicinity  of  the  first  harmonic  of  the  dominant  wave  frequency 
and  fall  into  agreement  only  for  frequencies  above  (approxi- 
mately) 7 cps. 

There  are  several  effects  brought  to  light  by  this  data. 

The  results  in  Figure  12  show  that  as  the  current  is  increased 
in  the  favorable  sense,  the  spectrum  energy  level  decreases 
significantly,  here  by  an  order  of  magnitude.  This  could  per- 
haps be  the  result  of  a reduced  shear  velocity  produced  by  a 
relaxing  surface  when  the  current  is  in  the  direction  of  the 
wind . The  adverse  current  does  not,  however,  produce  the  oppo- 
site effect.  The  principal  effect  of  an  adverse  current  has 
been  shown  in  Figure  20  to  be  an  increase  in  the  spectrum  energy 
only  in  the  vicinity  of  the  harmonic  of  the  dominant  wave. 

These  results  taken  together  suggest  the  effect  to  be  rather  of 
['etch  (or  time)  in  that  the  spectrum  has  more  time  to  develop 
in  the  case  of  an  adverse  current,  owing  to  the  correspondingly 
smaller  mean  convective  speeds.  Comparing  the  results  in  Fig- 
ure 13  for  the  flat  beach  to  those  in  Figure  19  for  the  inclined 
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beach,  at  the  same  wind  speed  and  very  nearly  the  same  local 
currents,  the  spectrum  energy  reduction  due  to  a favorable  cur- 
rent is  significantly  larger  for  the  inclined  than  for  the  flat 
beach.  This  suggests  one  possible  effect  of  local  current 
gradient . 

5.  CONCLUDING  REMARKS 

An  understanding  of  the  individual  and  combined  effects  of 
wind  speed,  current  speed  and  direction,  current  gradient  and 
fetch  (or  time)  in  win! -wave -current  interactions  is  of  obvious 
and  immediate  importance.  The  emphasis  in  the  first  experiments 
has  been  rather  diffused  over  all  of  the  parameters  in  order  to 
gain  some  insight  into  their  relative  importance  for  more  ef- 
fective subsequent  experimental  attacks.  In  spite  of  the 
hinderances  imposed  by  a mean  flow  suffering  from  unwanted  3-D 
distortions,  the  results  of  the  present  experiments  provide 
clear  evidence  of  the  character  of  the  interaction  phenomenon 
as  well  as  useful  information  on  the  trends  produced  by  varia- 
tions of  both  wind  speed  and  current.  The  important  effects, 
distortions  and  trends  observed  in  these  first  experiments  have 
been  described  and  briefly  discussed.  Much  more  data  is  needed, 
however,  before  a clear  picture  of  the  processes  involved  and 
the  roles  of  each  of  the  parameters  can  be  brought  to  the  front. 

We  expect  that  the  effects  of  wind  speed  and  current  will 

be  correlated  with  the  speed  ratio  /U  . To  minimize 

1 D0  c max 

mean  flow  distortions  due  to  sources  other  than  direct  wind- 

wave  -current  interactions,  the  values  of  and  U must  be 

5 Dq  c max 
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judiciously  chosen.  Recall  that  theoretical  considerations 
suggest  a speed  ratio  (absolute  value)  in  the  vicinity  of  unity 
for  maximum  interaction  effects.  Small  speeds,  on  the  one  hand, 
can  produce  conditions  where  the  effects  of  interactions  are 
also  small  and  difficult  to  perceive;  while  large  speeds,  on 
the  other,  can  produce  conditions  where  travel,  and  hence,  in- 
teraction times  are  short  in  a length-limited  facility.  Future 
experiments  performed  at  a moderate  wind  speed,  say,  4.75  m/sec, 
with  values  of  current  selected  for  a broad  range  of  speed 
ratio,  would  provide  data  to  elucidate  the  role  of  the  speed 
ratio . 

It  is,  no  doubt,  significant  that  over  the  range  of  condi- 
tions tested,  the  distortions  within  the  drift  layer,  supposedly 
produced  by  the  presence  of  the  current  field,  and  the  corres- 
ponding modulations  of  the  surface  wave  spectra  could  no  where 
be  called  other  than  very  small.  These  observations  will  be 
cast  into  more  quantitative  measure  when  the  RRI  results  become 
available.  Distortions  within  the  drift  layer  have  been  mea- 
sured. The  corresponding  temporal  frequency  spectra  results 
presented  here,  however,  show  no  striking  modulations,  either 
narrow  or  wide  band. 

Before  the  laboratory  results  can  be  applied  to  conditions 
encountered  in  the  field,  the  appropriate  scaling  laws  must 
first  be  determined.  How  to  compare  the  laboratory  drift  layer 
to  a typical  oceanic  drift  layer?  How  to  compare  the  effective 
internal  wave  period  and  wavelength  in  the  laboratory  to  those 
for  a typical  oceanic  internal  wave?  How  to  compare  the  measured 
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results  for  a stationary  beach  to  those  for  a progressive  in- 
ternal wave?  These  questions  are  presently  under  theoretical, 
consideration,  while  the  findings  and  results  of  the  experimental 
tests  are  being  carefully  analyzed  and  applied  to  the  design  of 
a modified  facility  and  the  next  test  program. 
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FIGURE  4 - WIND-INDUCED  AND  WAVE-INDUCED  SURFACE  CURRENTS 
UNDER  VARIOUS  WIND  VELOCITIES 
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FIGURE  6 - RELATIVE  DRi FT  VELOCITY  WITH  CURRENT  - WITHOUT  CURRENT 

VS. 

DOWN  STREAM  DISTANCE  FROM  LEADING  EDGE  OF  BEACH 

( From  Figure  5 ) 
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FIGURE  8 - DRIFT  CURRENT  VS  DEPTH,  WIND  SPEED  - 4.75  m/sec 
(Data  Set  B 02,  Flat  Beach) 
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FIGURE  10  - DRIFT  CURRENT  VS  DEPTH,  WIND  SPEED  - 7.65  cm/sec 
( DATA  SET  B02,  FLAT  BEACH  ) 
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FIGURE  14-  DRIFT  CURRENT  VS.  DEPTH  AT  STN  6,  WIND  SPEED  ^2.75  m/sec 
( DATA  SET  B4,  BEACH  ANGLE  - 2.65°  ) 
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FIGURE  15  - DRIFT  CURRENT  VS.  DEPTH  AT  STN  6,  WIND  SPEED  =4.75  m /sec 
( DATA  SET  B4,  BEACH  ANGLE  = 2.65°  ) 
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FIGURE  16 -DRIFT  CURRENT  VS.  DEPTH  AT  STN  6,  WIND  SPEED  - 6.67  m / sec 
( DATA  SET  B4,  BEACH  ANGLE  =2.65  ) 
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FIGURE  17  - DRIFT  CURRENT  VS.  DEPTH  AT  STN  6,  WIND  SPEED  = 7.65  m / sec 
( DATA  SET  B4,  BEACH  ANGLE  =2.65°  ) 
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x (M),  DISTANCE  FROM  LEADING  EDGE  OF  BEACH 


FIGURE  18  - CURRENT  ALONG  BEACH.  (Data  Set  B4,  Beach  at  2 .65°) 
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FIGURE  19  - WAVE  HEIGHT  - ENERGY  SPECTRA  AT  STATION  6 
( DATA  SET  B4,  BEACH  AT  2.65°  ) 


FIGURE  20  - WAVE  HEIGHT  - ENERGY  SPECTRA  AT  STATION  6.  (From 
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